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We suggest an approach to the problem of free electron spin evolution in a semiconductor with 
arbitrary anisotropy or quantum structure in a magnetic field. The developed approach utilizes 
quantum kinetic equations for average spin components. These equations represent the relaxation 
in terms of correlation functions for fluctuating effective fields responsible for spin relaxation. In 
a particular case when autocorrelation functions are dominant, the kinetic equations reduce to the 
Bloch equations. The developed formalism is applied to the problem of electron spin relaxation 
due to exchange scattering in a semimagnetic quantum well (QW) as well as to the spin relaxation 
in a QW due to Dyakonov-Perel mechanism. The results permit to separate the longitudinal Ti 
and transversal T-j relaxations times in a strong enough magnetic field and to trace the cases of 
undistinguished parameters Ti and Tb in zero and small magnetic fields. Some new predictions of 
the developed theory are discussed. 



Ctf 



-a 
c 
o 
o 



> 

00 
m 

o 



c3 



1 

C 

o 
o 



X 



I. INTRODUCTION 

The achievements in the investigations of ferromag- 
netism (FM) of diluted magnetic semiconductors (DMS) 
W and in spin injection technology with the possibil- 
ity to fabricate the spin-controlled devices have sparked 
renewed interest to studies of spin relaxation of electrons, 
holes and excitons (see J3J and references therein). Dur- 
ing the recent years, the optical picosecond technique ac- 
quire ability to observe the electron and hole spin kinetics 
in semiconductor quantum structures over wide range of 
temperatures and magnetic fields Q , || . In spite of in- 
tensive work the spin relaxation processes are still not 
fully understood in these structures. One of the reasons 
is the deficiencies in the present theories. If, for example, 
the theory of giant spin splitting of the electron states in 
DMS has a good base in terms of molecular exchange 
field, |q], spin relaxation theory has, by now, a huge gap 
between phenomenological and microscopic approaches 
to this problem. 

The phenomenological description of carrier or exciton 
spin relaxation in semiconductors is commonly started 
with Hamiltonian Hd without dissipation terms respon- 
sible for spin relaxation. For a particle with spin S 
(electron, exciton, etc.), the number of equations for av- 
erage spin values Aj =< S^SyS™ > (subscripts de- 
note Cartesian components; k,l,m satisfy the condition 
1 < k + I + m < 2S ; symbol / marks all indexes in the 
right hand part of Ai) is finite and is determined by the 
equations of motion ( K = 1) 
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(i) 



should include the relaxation part Rj to the Eq. (Q). 
Generally speaking, the relaxation of each Aj to thermal 
equilibrium value A® might depend on full set of spin 
averages. Therefore relaxation part reads 



Ri 



E 



(Aj-A°j) 



T i,J 



(2) 



Sometimes, the number of relaxation parameters tJ j can 
be reduced due to the presence of specific physical mech- 
anisms (see, for example, J7J). Nevertheless, it should 
be mentioned that a general procedure (for estimation of 
magnitudes of all relaxation parameters) that allows to 
reduce the numbers of relaxation terms (||) is absent. 

Other problem arises when one tries to relate the pa- 
rameters rfj with the spin-flip rate W (this is the aim of 
most microscopic calculations based on different known 
relaxation mechanisms). To clarify this problem we con- 
sider Bloch equations for average spin components S^ in 
a magnetic field directed along OZ axis 




Sx/T 2 

S Y /T 2 
(S z - S%)/T x 



(3) 



The system (nf) describes the evolution of spin system 
without relaxation. To take it into consideration we 



Eqs (||) looks like Eq. (|l|) , (||) , where S z is a thermal equi- 
librium value of Sz, w is Zeeman splitting in a magnetic 
field. Eqs (|h involve two relaxation times, the longitu- 
dinal T\ and transversal T 2 . In the framework of phe- 
nomenological Eqs (0) it is not apparent how the param- 
eters T\ and T 2 having different physical meaning are 
related to spin-flip probability W. Moreover, the trans- 
formation from the case of T± — T 2 at zero magnetic field 
(w = 0) to the case T\ ^ T 2 at u 7^ cannot be traced 
both with the help of Eqs (0) and by the microscopic 
calculations of spin- flip rate W. 



The present work is an attempt to "fill the gap" be- 
tween the phenomenological approach to spin relaxation 
and microscopic calculations of spin flip rates in the case 
of electron relaxation with spin 1/2 in a semiconductor. 
To do that, we derive the quantum kinetic equations for 
density matrix of spin 1/2. The equations have been de- 
rived in terms of correlation functions of dissipative sub- 
system. Despite the presence of some symmetry of corre- 
lation parameters, the kinetic equations can be only re- 
duced to the form that remains, nevertheless, more com- 
plicated than Bloch equations. In the following sections 
the developed theory is applied to electron exchange scat- 
tering on the magnetic ions in a semimagnetic quantum 
well (QW). Then, Dyakonov-Perel mechanism is consid- 
ered with respect to anisotropy of relaxation times in a 
QW and their reduction in a magnetic field. We discuss 
also the possible applications of the results obtained here. 



II. BASIC EQUATION FOR DENSITY MATRIX 

Let the electron spins represent a small part of total 
system so that all other variables is related to the bath 
being at thermal equilibrium. Corresponding Hamilto- 
nian reads 



H = H S 



H L + H, 



SL- 



(4) 



Here H$, Hl, and Hsl are, respectively, the Hamiltoni- 
ans of electron spin, dissipative subsystem (assumed to 
be at thermal equilibrium with a bath at a temperature 
T), and their interaction. According to projection oper- 
ator method (Ref. B) the density matrix of the system 
(0) is expressed in terms of main part and the rest: 



where 



p = fa + v, 



a = Ti L p; 

exp(-/3ff L ) 
1 Tr L cxp(-/lff L )' 
V = Pp=(l- fTr L )p, 



(5) 



^=-iH s a-C(a)-W(p(0)); 



with 



C(a)=Tr L V f S (t,t')Vfa(t')dt', 
Jo 

V(p(0))=Tr L VS(t,0)Pp(0), 

S (t, t') = exp | -* J {H s (r) + H L + PV) dr\ . (8) 

Here the calligraphic letters mean the Liouville oper- 
ators introduced over operator space as Tiscr = [Hs,cr], 
etc. The term V (p (0)) vanishes for the initial condition 
p(0)=fa(0). 

The expansion of Eq. (||) up to second order in the op- 
erator of interaction yields the kinetic equation for a: 

^M = -tHsa - Tr L j V (t, t) V (t, t') fa (t) dt', (9) 
where 

V (i, t') = exp l-i j H s (r) dr \ x 



xV(t)exp}ifHs(T)dr\, 

V (t) = exp {iH L t} V exp {-iH L t} 



(10) 
(11) 



We consider the reduced density matrix a(t) at the long 
enough time scale (as compared to times of attainment of 
a thermal equilibrium of dissipative subsystem). Thus, 
the upper limit of integration in Eq. (KJ) can be extended 
to infinity. Under this assumption we obtain Markovian 
equation for a(t), providing that spin relaxation times T 
are sufficiently longer than correlation times r of dissipa- 
tive system responsible for spin relaxation. 

To rewrite Eq.(H) in the matrix form, we should use 
some base. As the base, here we use the eigenvectors of 
spin operator Sz ■ For S = 1/2, the most general form 
of interaction is: 101 



We are looking for the kinetic equation for reduced 
density matrix a with the Hamiltonian 



where 



H = H S +H L + V, 



Hs = H° s + (H SL ) ; 
V = H SL - (H SL ) 



(6) 



(7) 



so that (V) = 0, where (...) = Tr L / 

It can be shown that exact form of equation for the a 

is 



Hsl 



SJ7. 



while electron spin-Hamiltonian can be written as 



Hs 



L0S9 



(12) 



(13) 



u> is the spin splitting in a magnetic field. Eq. ( |13[ ) can 
be considered also as a definition of OZ axis as quan- 
tization direction. Substitution of Eq.(|12|) to Eq.(0) (at 
t — » 00) , permits to express the collision integral in terms 
of the different Fourier components of correlation func- 
tions for Q a . Their number can be reduced with respect 
to the spectral properties of correlation functions and to 



the following identities (valid for any operators A and 
B): 



{A{r)B) L 



= e -0u 



\B(r)A)_ 



= e -P» 



{BA(r)) u 



(14) 



where the definitions of A(t) and B(t) are similar to Eq. 
(pd|); B = 1/T. The definition of Fourier transformation 
of correlation function is 



(A(r)B) u = 



2tt 



(A(t)B), 



(15) 



It is convenient to introduce the average values of spin 
components instead of density matrix a: 



(Sx) 
(Sy) 
(Sz) 



X = TraSx , 
Y = Tra S Y , 
Z = TraSz- 



(16) 



Straightforward calculations with the aid of Eqs( [l4|Jl5| ) 
yield following system of equations for components of (S): 



(17) 



where relaxation matrix T is determined via correlation 
functions 





iiiv (w) = (fy, (r) tt v ) u 



(18) 



Equation (p7|) describes the relaxation of spin 1/2 to 



its equilibrium state X (t) 



Xo,Y(t) 



Yq, and 



Z(t) 



Zq, with Xq = Yq = 0, and 



1 Buj 
Zn = — tanh — . 

2 2 



(19) 



7T I 


^ 7L + ™7 TO -™7zy 


-nj X z 


»l 


-nixy 7° z + nj xx 


-njy Z 


V -n-fxz -njyz 


nidxx +7to) 



The general form of relaxation matrix T is quite com- 
plex. Below, we consider two important specific cases 
permitting to simplify the components of correlation ma- 
trix j^ (to) Jig ); v, /i = X, Y, Z. First, we assume that 
the matrix ( |18[ ) is symmetric, 7 M „ (w) = 7^ (u>). Thus, 
the magnetic field influence on Hamiltonian of dissipa- 
tive subsystem is disregarded (Ref. ||l0| ). With respect 
to symmetric property of correlation function, one can 
obtain 



(20) 



n = n(w) = (1 + e^) /2, 7m „ = 7/lv («) , 7°, = 7m- (0) ■ 
Notice that Eqs (J17|) , (p0|) are, generally spiking, derived 
for the case of anisotropic medium and arbitrary spin 
splitting uj. 

In the case of zero spin splitting (uj — 0) and isotropic 
medium, the relaxation matrix T (EG) becomes symmet- 
ric. In this case, Eq.(pfl) are similar to those obtained 
in the work H in terms of fluctuating internal magnetic 
field n. 

In the case of strong enough spin splitting, one can 
expect that j a a(uj) << j zz (0). If we introduce Tf = 
7TH (jxx(uj) +7to(w)) and T 2 _ = 7T7 ZZ (0), and omit the 
off-diagonal components of T, the latter inequality per- 
mits to rewrite the Eq.(0 ) to the form of Bloch equa- 
tions (a). 

Second important case corresponds to antisymmetric 
coefficients 7^ (ui) = — 7„ M (uj) for v ^ /1. In this case we 
can obtain other form of relaxation matrix: 



Tzz + n lyy + ni lxy 







nijy Z 

7° 2 + n~f xx + nij xy nij zx 

n(-f XX + -fyy + lJ X y) 



(21) 



In the limit of zero magnetic field, the coefficients 
7/jw (uj — > 0) vanish if v ^ /i; thus the Eq.(pi|) transforms 
to matrix ( f20| ) with zero off-diagonal components. The 
spectral properties ( |14| ) with respect to definition of ti(uj) 
permits to observe that both relaxation matrix ( pl| ) and 
matrix ( |20| ) are even functions of the electron spin split- 
ting uj: T(uj) = T(-uj). 

The Eqs @ with definitions ©,@ are the main 
result of theoretical background of the work. In deriva- 
tion of kinetic equations ( |17| ) we assume neither any spe- 
cific properties of heat bath with Hamiltonian Hl nor 
any specific interaction mechanism H e L- We shall show 



how Eqs ( \L7\ ) can be used in the particular cases corre- 
sponding to ( p0| ) and (21). The short review of the re- 
lations between correlation functions (|15|) and two-time 
temperature Green's functions (with well-known calcula- 
tion technique) is presented in Appendix. 



III. EXCHANGE SCATTERING ON MAGNETIC 

IONS IN THE 2D QW 



Consider the problem of electron spin exchange scat- 
tering on the magnetic ions in a 2D quantum well (see 



Ref. |ll| jT^]). The carrier-ion spin-spin interaction op- 
erator has the form of following contact interaction 



H 



SL = -a^2 SJs e<5( r - R j 



(22) 



where a is exchange integral, S J is the magnetic ion spin 
at the site Rj; s e and r are the spin and coordinate 
of electron respectively. The sum is performed over all 
magnetic impurities. The Hamiltonian of dissipative sub- 
system include Zeeman energy of magnetic ions and the 
electron kinetic energy s^ with in-plane wave vector k. 
For simplicity sake, we restrict our consideration by the 
ground electronic confinement state only. So, the Hamil- 
tonian is assumed to take the form: 



Hr 






UqS z 



^ £ka Lr a k^ 



(23) 



where ojq is the magnetic ion Zeeman splitting in a mag- 
netic field directed along OZ axis; a^ and Ok, ff are the 
creation and annihilation operators; spin index a is in- 
troduced to normalize the chemical potential for total 
numbers of electrons. In our case 



H\ a U ak .<7 ='%2fk,* = 1 - 



(24) 



k.a 



k.a 



Now we can perform the renormalization procedure (m) 
with operator (23) and represent (second-quantized) in- 
teraction V (J12 ) in the form 



k^k' 
., j(k-k')R| t 



(25) 



ft; 



s. 



EE^i^r^ 2 ^-^^ 



k^k'.cr 

k^k' 



Here So is an area of a sample, i/)_l(Rj) is a perpendicu- 
lar (to a plane of a structure) component of confinement 
wave function. The term with k =/= k' is excluded from 
Eq. (Ea) due renormalization procedure (In). This proce- 
dure yields following Hamiltonian of dynamic subsystem: 



Hs 



ios ez = {(jj e + G eZ )s e 



(26) 



Eq.(p6[) defines spin splitting u in a total field. This 
field consists of external magnetic field (with Zeeman fre- 
quency uj e ) and exchange molecular field 

Ge = -|;E< sJ >W)r = 



So-T 
J 



a(S> / n m {z')\i>{z')Y dz' . 



(27) 



Here n m (z') is a local concentration of magnetic ions 
that is assumed to be the function of coordinate z di- 
rected along the growth axis Z' . In the case of magnetic 
QW with width L w and non-magnetic barriers, the inte- 
gration should be performed between —L w /2 and L w /2, 
while n m {z) should be replaced by average concentration 
n m in the QW. 

We are looking for the Fourier transformation of 
Green's function 

which (according to Eq. (|25|)) equals to 



7„., (") = §* E E Ei^fe)i 2 i^feoi 2 e< k - k O^+K kl -OH,^w ;kiCTi;k , CTiMj 

k,k',o-ki,k^,cri j,j' 



where 






^"^"k',^ (*) ; S£ <,.,«,! °ki,<ri 



(28) 



The explicit form for equation of motion (A4) for the cor- 
respondent Green's function wGjf'" (/a, v = X, Y, Z; q is 
a rest of indexes that determines the correlation function 
(pq)) is approximated by the following system of equa- 
tions 



^•,j'0k,ki0k',ki<Jo-,<T' /c . 



«cf - x («) = ""•"»■»■!;"• ^- g ( S %) [/ k)CT (i - / kV ) + / k ',.(i - / k ,.) 



2tt 
+ (ei* - £k) Gf •* («) - ^ G^ X (w) 



*<#*(«)=< 



^jj'^k.ki^k'.ki^Vcr' (5z) 



2vr 



[/k> - /k. CT ] + (£k< - £k) G*' A (w) + iw Gf A (, 



(29) 



(30) 



where we use (SxSy) = — (SySx) = i (Sz) /2; the oc- important for calculation of I\ So, taking into account 
cupation numbers /k,o- are determined by Eq.(E4f). Ac- that u> — ui + ie with infinitesimal e. one can use following 
cording to Eq.(|A3j), the imaginary part of the G%' v is identity in the solution of Eqs ( p9| , p0| ): 



u - (e' - e) 



[uj-(e>-e)Y-u>$ 



-i- [5 (w - (e' - e) - w ) + 6 (u - (e' - e) + w )] 



7 , , ° 2 o "^ ~* o ^ ^ ~ ( £ ' - £ ) - w o) - 5(w - (e' - e) + Wo)] 

[w-(e'-e)] -wtf ^ 

The final result for Fourier transformation of correlation function looks like 



(31) 



2n 



K x,r,x,j' , -, 

k 1 (7;k / ,cr;ki,cri ;k^,tTi \ / 
k^jk'jO-jkijO-ijkj.o-i ^ n 



(Gf*( w +*e)-Gf> x (w-*e)); 



K 



e>x («) 



^"^ k ' /k '- kl< ^{<^) [/k, CT (l - /kv) + /k>(l - /k,.) 



x [5 (w - (e' - e) - w ) + 5 (w - (e' - e) + w )] + 

+ ^ [/k> - / k . ff ] [^ (w - (e' - e) - «o) - J (« - (e' - e) + w )] 
2n 



(32) 



lf £;fcki l oi;ki 1 oi M = 4^^,j'^k J k' 1 ^k',k 1 ^, - 1 j —£~ [h',a - fk,cr] * 

x[5(lu- (e' - e ) - Wo ) + 5 (w - (e' - e) + w )] + 

+ ^[/k, ff (l - h>,a) + /k', ff (l - h,a)] [S (w - (e' - e) - wo) - ^ (w - (e' - e) + wo)]}. 



(33) 



with e = £k,e' = £k' 



We assume that neither 2D structure stress, nor 
internal magnetic field affects the transversal mag- 
netization (which means \S X ) = (S Y )). Thus 



k : CT;k : CT;ki ,cji ;k^,cri V / 
-X,j;Y,j' 



K, ,3 \ ,' 3 , , , (uj) and 

'Y,j;X,j' 



K'^Cr u v M = -if' ™ k t , (w). Lat- 
ter equality shows that the exchange scattering corre- 
sponds just to the case of relaxation matrix (pl|). Direct 
calculation of off-diagonal matrix elements shows that 
K*> j 'Z' 3 \ v M = KY' j 'lf\ ,, («) = 

k,CT;k ,CT;ki,fTi ;k^,(Ti v / k,a";k .<7;ki,(7i ;k^,(Ti v ' 

since (Sx) = {Sy} = 0. In a similar manner, one can 
find 



-Z,j;Z,j' 



1 



K Z'X-M^;K^ (0) = 2 S ^' 6k ^ 6k '^ 6 ^ ( S ^ X 
x [/ k , CT (l - /k>) + /k'..(l - /k, CT )]<5 (e' - e) . (34) 

The relaxation parameters in the matri x (|2l| ) can now 
be obtained with the help of (p3,p3). Eq.([24|) reads now 



X^ k > CT 



2irh 2 
SnmT 



,-Ek/T 



(35) 



where m is the in-plane effective mass of the carriers. 
Substitution of Eq.(p4|) to the formula irn {Viz (t) Qz) 
gives following result in the case of /k,<j <C 1 : 



(k#k') 



7i"7° z 



*§? E Y.( S z)\^)\ 4 lfK« + h'A6(e-e') 



k,k',cr j 



Calculation of sums in Eq.(p6|) yields 

(Sl), 



IT n a 2 m 

h Jzz 2L w h? 



(36) 



(37) 



where dimensionless parameter / reflects the overlap of 
^-function and magnetic ions: 



I = L„ 



"^l^z'f dz> 



(38) 



In the case of infinitely deep QW containing the mag- 
netic ions with concentration n m , the integration in pa ) 
gives I = 3/2. 



To obtain other elements of the matrix (21) we should 
find j xx and ^ xy according to Eqs(p2|),(|33|). It can be 
shown that contribution of a second term in Eq.(32) is 
less then that of first one, while j xy ^C j xx — j yy , so that 
we can find 



(k#k') 



kti~{ xx 



-*Sr E E 



(S 2 x) 



l^(^)| 4 [/k, CT + /kVp(w - (e - e') - uj a ) + S(lu - (e - e') + w )]. 



(39) 



k,k',cr J 



Sums in Eg. (J39|) are evaluated similar to those in 
Eq.@: 



n 



o?m 



(s%) 



F 



/W Wo 

\T' T 



(40) 



L^ 3 ™ 2 
where we introduce the function 

Ffoy) = - [2 + exp (-|a; + y\) + cxp (-|.t - y\)\ (41) 

that decreases monotonically from 1 to 1/2 and reflects 
the suppression of electron spin fluctuations by a mag- 
netic field. Thus, the kinetic equations (plf) transforms to 
the Bloch equations (|3j) with following relaxation param- 
eters 



with the help of Eqs @, @, @, @. One can see 
that magnetic field suppresses the longitudinal relaxation 
while it enhances transversal relaxation for S > 1/2. In 
the limit of saturated magnetic field, uj,uj » T, the de- 



creasing of Ti reaches the factor 



2 5+1 



Note, that the 



3 S+l/2 - 

decreasing of transversal relaxation time in a magnetic 
field was observed for Mn spins but was not explained 
till now. 

At the end of this section we estimate the relaxation 
time of the electron in C d\- x Mn x T e QW with x=0.017 
and L w — 8 nm. At zero magnetic field one can find 
t p = 5 ps that is quite close to experimental value ob- 
served recently via time resolved magneto-optic Kerr ef- 
fect. [HI 



Ti J r L w h 3 '' 



{Sl)F(^ 



LU LdQ 



T" T 



(42) 



IV. DYAKONOV-PEREL MECHANISM 




(43) 



The right hand parts of Eqs fl37|) and ([40|) equal to those 
of Eqs ([42| ) and (J43|) respectively at zero magnetic field as 
it should be expected. Moreover, as w — > 0, ojq — > 0, Eqs 
([42]) , (|43|) transform just to the double rate W of spin flip 
relaxation obtained with Golden Fermi Rule, |lll|| , |12| 

With magnetic field increasing, one can obtain in the 
limit w,wq >T, 



Ti 


1 




= 25+- 


T 2 


2 



(44) 



Eq. ( p4[ ) shows that rate of phase relaxation exceeds the 
longitudinal relaxation rate with factor about one order 
of magnitude in the case of ions Mn 2+ (S — 5/2). In 
the case of equidistant Zeeman splitting uiq of the mag- 
netic ions, the magnetic field dependence of spin square 
averages is: 



(^)=icoth^b s (^), 
<S|) =S(S+1)-2(S 2 X ), 



(45) 



where bs(x) = — {Sz) = SBs(Sx) is non-normalized 
Brillouin function: 

b s (x) =(S+l) coth(5 + \)x-\ coth I (46) 

Figure 1 shows the magnetic field dependence of the 
longitudinal and transversal relaxation times calculated 



In a number of works |Q, 15 1, |[|, it have been shown 
that a fluctuating effective magnetic field accompaning 
the electron momentum scattering in the semiconduc- 
tors with non-centrosymmetric lattice results in rather 
efficient spin relaxation. In this section we show how 
Dyakonov-Perel mechanism can be incorporated into for- 
malism developed above. 

Hamiltonian of dissipative subsystem is assumed to 
have the form 



H L =Yj £ k«k,, 



flk., 



k,a 



E 

k,k',<r 



U 



k,k' a k,o- ak ' ! °'- 



(47) 



Second term in Eq.(|47J) has symbolic meaning: we as- 
sume Vk,k' to be an operator responsible for electron scat- 
tering by the impurities, phonons, or electron-electron 
collisions. We should also introduce the electron scat- 
tering changing randomly an effective magnetic field 
B = b(k)/g e /x that corresponds to band spin splitting 



Hi. 



^ s e b(k)aj c a k> , 



k,« 



Therefore, the operator of this field reads 
tin = J2b fi (k)al a a Ka . 



(48) 



(49) 



The specific form of the function b^ (k) is determined by 
a band Hamiltonian (see, for instance, Ref. M) and de- 
pends on dimensionality of the system under considera- 
tion (Ref. H|). 



We introduce the Green's function 



G = ^2,3 — 



a k a ak »! a ki° k i 



(50) 



For simplicity, here we suppress spin indices while 1 



ki, 2 = k 2 , etc. Equation (A4) with Hamiltonian ( |47| ) 
assumes the form 



[w - (£i - e 2 )\ G 2 , 3 = ^M/ X (l - h) + 



/ y ( V3 ! 4G , 2,4 — V^Gi^ 



(51) 



We are looking for Gi : i that can be represented via 62,3. 
Solution of the Eq.|3l| depends strongly on k, k' - depen- 
dences of operators Vk : k'- We assume that following gen- 
eral relation for Green's function holds (Ref. |l8||) 



G1.1 — 



A(l-/i) 
■k (w + iTi t i) ' 



(52) 



where Tu means the rate of electron moment relaxation 
due to scattering processes described by second term in 
the Eq.(|47|). Thus, the relaxation parameters are defined 
by 



n ni ^ v 



r k , k /s 



L^w^w ^^ 



h 2 o 



(53) 



•k,k 



Let us consider the electron spin splitting caused by 
k 3 - terms in band Hamiltonian (Ref. Urn). This mecha- 
nism had been considered for 2D-electron gas in the work 
Jl5| . If external magnetic field is directed along the axis 
of structure growth, one can obtain 



w 



with 



v/r ! 



m 3 / 2 (2£ 9 )V2' 



K x — fcxQ ) ^y — fcyQ 5 K z — U, 



(54) 



(55) 



(56) 



where m is in-plane effective mass, q 2 is defined via elec- 
tron ip- function confined in a QW, 



^ 



dz 



V- 



(57) 



We see that only two components of relaxation param- 
eters, 7^2 and j Vj y have nonzero values. The calculations 
with the help of Eq. (B3) give rise to following expression 
for transversal relaxation time 



1 



-n^ XtX 



-n"fy, : 



a 2 K z q A T 



2m 2 E g 1 + (wr e ) 2 



(58) 



where the electron scattering time is defined via averaged 
value of r k k : 



r- 1 = h- 1 (r k , k ) 



(59) 



Note that approximation (|59j) is not nessessary for cal- 
culations of T\ and T%. The dependence r e (k) can be 
taken into account for more accurate theories (see Ref. 
§)• 

According to Eq.(|2(j), the longitudinal relaxation rate 
is 



'-- n[pj x ^ x + "/y,y ) rr, ■ 
1 1 n ±2 



(60) 



We can see that Tican be shorter than T^ in the systems 
with strong anisotropy . 

Other important case corresponds to magnetic field di- 
rected along OZ axis in QW plane while OX axis is a 
perpendicular to this plane. One can find 



V] Ky 



-k y q 



k z q 2 . 



(61) 



In this case, three parameters describe electron spin 
relaxation: 



1 
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a 2 K 2 q 4 T^ 
2m 2 E g 

aWq*T r 
2m 2 E g ' 

a 2 H 2 q 4 T 



1 + {uJT e f 



2m 2 E g 1 + (ujT e y 



(62) 



We can see that transversal relaxation (Eq.[58|) can be 
suppressed significantly by the strong enough magnetic 
field, Lor e >> 1, due to dynamical averaging effect in the 
case of Faraday geometry of experiment |l6[ . In the case 
of Voigt configuration (Eqs (p2|)), magnetic field can sup- 
press the transversal relaxation time by no more than two 
times. The effect of drastic magnetic anisotropy of spin 
relaxation in QWs can be used to recognize the specific 
mechanism under consideration. Actually, the typical 
scattering time of an electron in semiconductors is about 
r e ~ 10~ 12 . To make the magnetic field anisotropy ef- 
fect visible, one should apply quite strong magnetic field 
to non-magnetic structure. On the other hand, consid- 
ered anisotropy effect can be amplified in DMS structures 
giving rise to inequality u> » luq at moderate magnetic 
fields (Eqs ©, @). 



V. CONCLUSIONS 

In this work, we develop the microscopic theory of elec- 
tron spin evolution in semiconductors in terms of quan- 
tum kinetic equations for arbitrary mechanism of spin 
relaxation. The relaxation term in these equations has 
been derived as a matrix of correlation functions provid- 
ing relatively short correlation times ( as compared to 
spin relaxation times, r << T) of dissipative subsystem. 
Our theory permits to distinguish the phase and energy 
relaxation processes and to capture the case of zero and 
small magnetic fields. 



relaxation in QW in terms of Eqs fll7| , |20| , |21| ) promises 
new insight into this problem. The investigations of a 
hole relaxation in QWs or uniaxial semiconductors will 
be reported elsewhere. 
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APPENDIX A: 

Equation u\j\) reduces the problem of spin relaxation 
to calculations of correlation functions of operators J~2 re- 
lated to a heat bath at a thermal equilibrium. Theory of 
correlation functions for quantum-mechanical operators 



2 3 4 

B, Tesla 



is well developed (see e.g. Ref. 18 1). In this Appendix, 
we present few equations that are needed in subsequent 
analysis of relaxation parameters. 

Introduce first two-time temperature retarded Green's 
function (Ref. @) 

G(t, f) = «n„; n„» = -ie(t - o ([n p (t) , ti„(t)]) . 

(Al) 

The bracket [fip (t) , ft v (t)] denotes a commutator or an- 
ticommutator for operators f^ (t) and Q v (t). Fourier 
transformation of Green's function (Al) reads 



FIG. 1. Magnetic field dependencies of the rates of longitu- 
dinal (i — 1) and transversal (i — 2) electron spin relaxation 
at T = 4K in Cd 1 ^ x Mn x Te (x = 0.017) crystal. 

We show that spin relaxation process due to exchange 
scattering in a semimagnetic QW determines the longitu- 
dinal and transversal spin relaxation times. Specific cal- 
culations of relaxation parameters show that only longi- 
tudinal relaxation can be described by flip-flop processes 
which is suppressed in a magnetic field. On the other 
hand, transversal relaxation is due to the effective ex- 
change field fluctuations, which increase with magnetic 
field. Qualitatively, the latter effect has the same reason 
as growth of line width of spin flip Raman scattering by 
shallow donors in diluted magnetic semiconductors fll7| ]. 

The theory has been successfully applied to Dyakonov- 
Perel mechanism of carrier spin relaxation in a QW. The 
significant difference in transversal relaxation rates for 
parallel and perpendicular (to growth axis) magnetic field 
orientations has been found. 

Proposed approach to relaxation can be extended to 
any two-level system described by fictitious spin 1/2. 
From this standpoint, the analysis of heavy hole spin 



GH = ((^;^)) 



1 

2^ 



G(t, 0) exp(iut)dt. (A2) 



Fourier transformation of correlation functions (Eq.(15)) 
can be expressed in terms of that for Green's functions 
Eq.(P): 



. . .,. G(uj + ie) -G(u-ie) 

7Mw) = ihm —— . 

£->o 2n{u!) 



(A3) 



So, the problem is red uced to explicit calculation of 
Green's functions in Eq.(Al). The equation of motion 
for their Fourier component is |lq]: 



wG(w) = i- <[«„, OJ) + (((SVfli - H L %) ■ Q u )) u 



(A4) 



The solution of equation (A4) permit to find the correla- 
tion function, which, in turn, determines the spin relax- 
ation according to Eqs ([i"5|), (|2C|) and (pil"|). 
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